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Abstract 

Wind tunnel measurements of drag force at low 
speeds were taken for a laminar-flow body of revo- 
lution free of support interference. This body was 
tested at zero incidence in the 13-inch Magnetic Sus- 
pension and Balance System (MSBS) at the Lang- 
ley Research Center. The primary objective of these 
tests was to substantiate the capabilities of the MSBS 
for measuring drag force. The drag-force calibra- 
tions and wind-on repeatability data provided the 
means of assessing these capabilities. Additional in- 
vestigations included (1) the effects of fixing transi- 
tion, (2) the effects of fins installed in the tail, and 
(3) surface flow visualizations using both liquid crys- 
tals and oil flow. Two simple drag prediction codes 
were used to assess their usefulness in estimating 
overall body drag. These theoretical results did 
not compare well with the measured values because 
of incorrect or nonexistent modeling of a laminar 
separation bubble on the body. 

Introduction 

Support interference is a serious problem when 
testing models in wind tunnels (refs. 1 and 2). Fig- 
ure 1 shows three examples of model support inter- 
ference problems. In each case, the sting support re- 
quired a model geometry different from the aircraft 
geometry. The presence of the sting also distorted the 
base flow in each case. This combination of geome- 
try modification and distorted base flow produced 
the errors shown. Using a magnetic suspension and 
balance system (MSBS) is the only way to eliminate 
support interference completely. With the elimina- 
tion of the support, not only will the flow distortion 
produced by the sting be eliminated, but also many 
other advantages will accrue, such as 

1. elimination of model modifications to accom- 
modate the sting 

2. ease of model movement for dynamic testing 

3. fast, efficient testing at any attitude 

4. improvement in productivity by elimination of 
stings and struts 

Other potential advantages are presented in 
reference 3. 

An MSBS is a device capable of both suspend- 
ing a model in space and measuring the forces and 
moments acting upon it. These actions are accom- 
plished through the use of controlled magnetic fields 
interacting with a magnetized core placed in the 
model. In a wind tunnel, an MSBS uses electromag- 
nets positioned external to the test section. The core 
is either a permanent magnet, magnetized soft iron, 
or a solenoid within the model. The position of the 


suspended model is inherently unstable, and a closed- 
loop feedback control system which incorporates a 
method of sensing the model position and attitude is 
therefore required to stabilize the model by actively 
controlling the applied magnetic fields (refs. 4 and 5). 
The MSBS can counteract a range of external forces 
and moments applied to the model while maintain- 
ing the model at any arbitrary position. The amount 
of electrical current flowing in each of the external 
electromagnets is used to determine the forces and 
moments acting on the suspended model. 

The French at ONERA first demonstrated a wind 
tunnel MSBS in the mid-1950’s (ref. 6). Since 
that time, considerable research has been devoted 
to MSBS’s and many improvements have been made 
(ref. 7). Today, applying MSBS’s to routine wind 
tunnel testing of aircraft models is practical; how- 
ever, they are still not widely used. This is princi- 
pally because of the small size of all existing systems. 

Recent technical advances have removed some of 
the barriers to building large MSBS’s. One example 
is progress with large superconducting electromag- 
nets which would be needed to support the weight 
of the model and also to counteract the large aero- 
dynamic loads. Another advance is the use of a 
superconducting solenoid as a model core (ref. 8) to 
produce a model with the highest possible magnetic 
moment. Using this type of core would reduce the 
overall system cost by minimizing the required size 
of the external electromagnets. In addition to these 
advances, feasibility and conceptual design studies 
(refs. 9 to 11) have produced innovative MSBS de- 
signs. These new designs significantly reduce the 
cost of large systems, making them more affordable 
for large wind tunnels. Additional developments are 
needed in a few hardware technologies in order to 
fully use a large MSBS. These technologies include 
(1) very-high-capacity, energy-efficient power ampli- 
fiers and (2) wide-ranging, accurate, and versatile 
model position sensors. It is anticipated that work in 
these and related branches of engineering will result 
in steady development. 

Since the late 1970’s, in-house activity with 
MSBS’s has steadily increased at NASA Langley Re- 
search Center. In 1979 the U.S. Air Force loaned 
the 13-in. MSBS at the Arnold Engineering Develop- 
ment Center (AEDC, ref. 12) to researchers at Lang- 
ley. The system was initially used to gain operational 
experience with a working MSBS. Ownership of the 
system has since transferred to NASA from the Air 
Force. 

Early in 1984 this 13-in. MSBS was combined 
with a small, low-speed (M < 0.5) wind tunnel. The 
primary goal was to gain some practical experience 
with this system, with the emphasis placed on using 


and understanding the MSBS rather than on obtain- 
ing aerodynamic data. The 13-in. MSBS was then 
modified in several respects. The original AEDC 
X-ray position sensors were replaced with an electro- 
optical system which uses solid-state linear photo- 
diode arrays, and the original analog control system 
was replaced with a digital controller which uses a 
minicomputer (ref. 13). These and other improve- 
ments make the system easier to use. 

However, several limitations still exist with the 
13-in. MSBS. One is a limit on available angles-of- 
attack (±5°). Another is difficulty in suspending 
irregular-shaped three-dimensional models. Both of 
these limitations arise because of the design and 
nature of the electro-optical position sensing system. 
Despite these limitations, the emphasis at the present 
time is on using the 13-in. MSBS to obtain useful 
aerodynamic data. The existing position sensing 
system works well with most axisymmetric shapes. 
Therefore, a program is underway to obtain accurate 
drag data on various axisymmetric shapes at zero 
incidence. 

As part of this program, wind tunnel measure- 
ments of drag force at low speeds were taken on a 
laminar-flow body of revolution free of support in- 
terference in the Langley 13-in. MSBS. The primary 
objective of these tests was to substantiate the drag- 
force measuring capabilities of the MSBS. Generally, 
a multicomponent strain-gage balance is accurate to 
±0.5 percent of full scale. The basic accuracy of a 
multicomponent magnetic balance should equal or 
exceed that of a standard multicomponent strain- 
gage balance; however, the obvious advantage of an 
MSBS is the total elimination of support interference. 
In this context, aerodynamic data from an MSBS is 
more accurate than data obtained from a conven- 
tionally supported model with a strain-gage balance. 
The axisymmetric body used in this study was chosen 
because the unique data set produced would be of in- 
terest to researchers involved in designing minimum- 
drag fuselages. 

Symbols 

A polynomial coefficient from 

drag-force calibration equation 

A max maximum cross-sectional area 

of model, in 2 

b tunnel width, in. 

B polynomial coefficient from 

drag-force calibration equation 

C polynomial coefficient from 

drag-force calibration equation 


Cd 

drag coefficient 

Cd,o 

drag coefficient at zero lift 

Cd, o 

afterbody drag coefficient at 
zero angle of attack 

Cf 

skin-friction coefficient 

Cm 

pitching-moment coefficient 

Cn T 

yaw damping coefficient 

D 

drag force, gf 

h 

tunnel height, in. 

k 

incidence angle of tail, deg 

I 

electromagnet current, A 

L 

model length, in. 

l 

axial length, in. 

M 

free-stream Mach number 

P 

pressure, lbf/ft 2 

R 

Reynolds number 

r 

body ordinate, in. 

S 

reference area, in 2 

s 

distance along body contour, 
in. 

T 

temperature, °R 

t 

maximum model thickness, in. 

V 

free-stream velocity, ft/sec 

V 

local surface velocity, ft /sec 

X 

body abscissa, in. 

a 

angle of attack, deg 

0 

= y/\ - W 

A 

change in value 

e 

blockage factor 

A 

body shape factor 

T 

tunnel shape parameter 

Subscripts: 


B 

due to buoyancy 

b 

due to blockage 

c 

corrected value 

cv 

measured in the cavity at the 
rear of the model 

L 

based on model length 

0 

based on zero axial load 
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s static condition 

s' due to solid blockage 

t total condition 

w due to wake blockage 

wet based on wetted surface area 

wo wind on 

woff wind off 

Description of 13-in. MSBS 

Figure 2(a) shows a sketch of the arrangement 
of the 13-in. MSBS which was originally constructed 
at the U.S. Air Force AEDC in the early 1960’s. 
For this system, motion is controlled in 5 degrees of 
freedom with no provision for generation of controlled 
magnetic roll torque on the model. Figure 2(b) 
shows a photograph of the system. This system has 
four electromagnets, arranged in a V-configuration 
above the wind tunnel test section, that provide the 
lift force, pitching moment, side force, and yawing 
moment. The 13-in. MSBS has a lift force capability 
of about 6 lbf depending on the size and shape of 
the iron core in the model. The drag electromagnet 
counteracts the drag force. The test section for the 
tunnel passes through the drag electromagnet. 

All the electromagnets have conventional copper 
windings. The drag electromagnet is fed from a bipo- 
lar thyristor power supply, while the other electro- 
magnets are fed from a mix of unipolar thyratron 
supplies and rectified motor-driven variable trans- 
formers. The bipolar supply used with the drag elec- 
tromagnet enables it to both push and pull on the 
model. Note that each of the four main electromag- 
nets is composed of two sets of windings per pole, a 
set of bias windings and a set of control windings. 
The bias windings supply a steady current necessary 
to suspend a model. This reduces the amount of con- 
trol current needed during actual operation. 

The control system is completely digitized, with 
a minicomputer handling all the control functions 
(ref. 13). Typical feedback loop repetition rates for 
this system are on the order of 256 cycles per sec- 
ond. The 13-in. MSBS uses an optical position sens- 
ing system based on solid-state linear photodiode ar- 
rays. Laser light sheets are directed across the test 
section to illuminate the arrays. The model posi- 
tion and attitude are inferred by locating the shadow 
of the model on the arrays. This technique is very 
dependent on geometry. Abrupt changes in shadow 
positions are generally not tolerated by the control 
system; therefore, it is difficult to suspend many non- 
axisymmetric shapes, especially bodies with wings at 


various pitch or yaw angles. Reference 14 gives more 
details on this optical position sensing system. 

Figures 3, 4, and 5 are photographs of the 
laminar-flow body suspended in the 13-in. MSBS. 
Figure 3 is a view of the 13-in. MSBS from the con- 
trol room. Figure 4 is a close-up of the suspended 
laminar-flow body, and figure 5 is a view from the 
contraction section of the tunnel looking at the nose 
of the suspended body. 

Description of Wind Tunnel 

A small, low-speed (M < 0.5) wind tunnel was 
adapted to use with the 13-in. MSBS. As shown in 
figure 6, the tunnel is an open-circuit design with 
air drawn into the tunnel from outside the build- 
ing through a large bellmouth protected from out- 
side contaminants by a screen enclosure. Constant- 
diameter ducting extends the circuit to the first turn. 
Following the first turn, a quick (rapid area) diffuser 
and settling chamber with screens and honeycomb 
prepare the flow for entry into the contraction and 
test sections. The contraction ratio is 8.43. 

The clear General Electric Lexan thermoplastic 
test section is a modified octagon shape with major 
and minor axes of 12.56 and 10.69 in., respectively. 
(See fig. 7.) At the end of the diffuser, a section of 
2.5-in.-thick honeycomb is used to protect the fan. 
A turn identical to the inlet turn directs the flow 
to the fan section and then to the exit outside the 
building. A water-cooled, 200-hp, 6000-rpm, variable 
frequency electric motor drives the fan. Reference 15 
gives full details on the tunnel layout and calibration. 

The flow quality in the tunnel is as follows 

1. a maximum deviation of ±0.25 percent in dy- 
namic pressure across the test section 

2. velocity fluctuations in the streamwise direc- 
tion of 0.2 percent ±0.1 percent for the 
Reynolds number range of 2 x 10 5 ft -1 to 
2 x 10 6 ft -1 

3. flow angularity in pitch of about -0.5° relative 
to the cathetometer reference (see Test Proce- 
dure section) 

Description of Model 

Figure 8 shows a photograph of the model used in 
this study which has a long favorable pressure gradi- 
ent forebody. The body profile was generated from 
an eight-parameter class of rounded-nose, tail-boom 
bodies described in reference 16. This class of bodies 
was developed to verify a method of shaping axisym- 
metric bodies to produce minimum drag in incom- 
pressible, nonseparating flow at zero incidence. An 
optimization scheme constrained to exclude turbu- 
lent separation is used to minimize Cp. This scheme, 
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which achieves drag reduction by body-shape modi- 
fication, results in bodies with extended runs of lam- 
inar flow. 

A sting-mounted version of this shape was previ- 
ously tested in a towing tank at the United States 
Naval Academy (refs. 17 and 18). The model tested 
by the Navy was slightly modified from the design 
shape in the tail region in order to accommodate the 
sting. (See fig. 9.) One of the goals of these tests was 
to quantify the hydrodynamic performance of this 
body. Additional computational results on this body 
are found in references 19 and 20. 

The model used for the present study was 12 in. 
long with a maximum diameter of 2.67 in. The 
model, which weighed about 2 lb, was made in two 
longitudinal halves with a fiberglass skin that was 
foam filled around an embedded soft iron core. This 
ultralow-carbon iron core was magnetized by the ap- 
plied magnetic fields. The core was 1 in. in diameter 
by 6 in. long and it was centered in the model. Be- 
cause of the lack of roll control for the 13-in. MSBS, 
the model was weighted to provide a constant roll 
position. This was done during fabrication by plac- 
ing a 20-g brass weight against the inside surface of 
the model’s skin. 

The eight nondimensional parameters listed be- 
low were used to generate this body of revolution 
according to the method described in reference 16. 


r n = 0.50 

Xi = 0.80 

11 

bi 

r; = 0.50 

x m = 0.5555 

Si = 2.0 

ki = 1.50 

t = 0.3 


The symbols used are defined in reference 16. 
Table I lists the coordinates of the model. 

Description of Surface Flow Visualization 
Techniques 

Liquid Crystals 

A liquid crystal surface flow visualization tech- 
nique was used to determine transition and separa- 
tion locations. This technique is discussed in detail 
in references 21 to 23. Liquid crystals have properties 
of both liquid-phase and solid-phase materials. They 
can display the optical properties of solid crystals by 
scattering light very selectively. These optical char- 
acteristics can be displayed by altering the molecular 
structure of the liquid crystals. The molecular struc- 
ture is primarily influenced by temperature and shear 
stress, which cause the liquid crystals to selectively 
reflect various colors. 


Liquid crystals with a wide temperature band- 
width were chosen for this study so that shear stress 
would dominate in determining the liquid crystal 
color response. This shear stress response can be 
used to visualize various boundary layer flows, in- 
cluding transition and separation. Liquid crystals 
indicate transition by changing colors in response to 
the large shear stress change through the transitional 
region. (See fig. 10, taken from ref. 23.) Separation 
is indicated by a change from “color” to “no color” at 
a point on a body because C j is zero in a separated 
region. 

For the present study, the models were painted 
flat black. As noted in reference 23, this provides the 
best background surface for visualizing liquid crys- 
tals. The liquid crystals were thinned with a solvent 
and applied by spraying a thin coat on the model sur- 
face. A good application results in a coating thick- 
ness of approximately 0.0005 to 0.0010 in. One coat 
could be used for several tunnel runs at different 
Reynolds numbers. 

Pigmented Oil 

A pigmented oil surface flow visualization tech- 
nique was also used to determine separation loca- 
tions. This is a type of oil flow which uses a mixture 
of titanium dioxide (Ti0 2 ) suspended in refined min- 
eral oil. A small amount of oleic acid is added as an 
anticoagulant. After the mixture was brushed on the 
black model surface, the model was placed in the tun- 
nel and suspended. The tunnel was then brought to 
test conditions, and after the excess oil was wiped 
off, the Ti0 2 was left deposited on the surface. 

Description of Drag Prediction Codes 

SANDRAG Code 

A simple and efficient drag prediction code, 
SANDRAG, was used to obtain a theoretical esti- 
mate of C D ' Wet for this body (refs. 24 and 25). SAN- 
DRAG calculates the flow field and drag of bodies 
of revolution at zero incidence in incompressible flow 
by using a combination of a potential-flow method 
and boundary-layer techniques. This code has been 
used successfully to predict drag for some simple 
axisymmetric projectile shapes. The goal for this 
study was to verify the usefulness of this code in pre- 
dicting drag for the more complicated axisymmetric 
shape used. An axial distribution of source and sink 
elements is used to form the body for the potential 
solution, with the equations applying to both open 
and closed bodies. The body used for this study is 
considered an open body because the surface does 
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not return to the axis at the tail. The laminar- and 
turbulent-boundary- layer solutions use momentum 
integral techniques which account for the effects of 
surface roughness. Transition is modeled by taking 
into account the effects of pressure gradient and 
surface roughness. Separation is also modeled by 
predicting separation locations and estimating the 
drag based on the pressure in the separated region 
(separation drag). In the laminar boundary layer, 
separation is based only on the pressure gradient, 
while in the turbulent boundary layer, separation 
is based on the value of the shape factor. Finally, 
the base pressure coefficient is estimated with data 
correlations from various afterbody shapes. 

Nakayama and Patel Code 

A simple code was developed by Nakayama and 
Patel to calculate the viscous resistance, or skin- 
friction drag, of a streamlined body of revolution 
at zero incidence in incompressible flow (ref. 26). 
This code calculates the flow field and drag by us- 
ing a combination of a potential-flow method and 
momentum-integral boundary-layer techniques. The 
Nakayama and Patel code has been used to calcu- 
late fairly accurate drag coefficients for simple airship 
shapes in cases where the transition point is known 
or is predicted accurately (ref. 26). Once again, the 
goal in this study was to verify the usefulness of 
this code in predicting drag for the more complicated 
axisymmetric shape used. 

Several empirical correlations are used to predict 
transition. In addition, the experimental transition 
location can be used. This code was primarily de- 
veloped to evaluate the flow in the tail region. The 
turbulent-boundary-layer calculations are made with 
Patel’s integral method (ref. 27), which accounts for 
thick turbulent boundary layers. Separation is not 
modeled by this code; in addition, this code does not 
provide estimates of base pressure drag. 

Test Procedure 

For this study, the model was tested at zero in- 
cidence. The pitch attitude, vertical position, and 
horizontal position of the model were verified with 
an optical cathetometer. This is a remote fixed de- 
vice used to measure relative vertical and horizontal 
displacements. (See fig. 11.) The pitch attitude, ver- 
tical position, and horizontal position were accurate 
to approximately ±0.02°, ±0.002 in., and ±0.002 in., 
respectively. The model pitch attitude was not cor- 
rected for the flow angularity mentioned earlier. The 
yaw attitude and lateral model position relative to 
the test section sidewall were verified by shining a 
laser light sheet longitudinally across the top of the 


model. It is estimated that the yaw attitude and lat- 
eral position were accurate to ±0.6° and ±0.1 in., 
respectively. 

The initial phase of testing consisted of a drag- 
force calibration. The calibration involved suspend- 
ing the model with no wind and loading it in the 
axial direction. For this type of calibration the elec- 
tromagnets above the test section suspend the model 
by supporting its weight. The drag electromagnet 
works to counteract the applied load. The loading 
was performed with the traditional string, pulley, and 
weight pan combination shown in figure 12. 

A data acquisition unit (DAU) and a microcom- 
puter were used to read and record the data. (See 
fig. 13.) The DAU was interfaced directly with 
the shunts which provide a means of measuring the 
current in the electromagnets. For a given weight 
or load, shunt voltages and magnet currents were 
recorded with the microcomputer. Voltages were 
converted on-line to obtain equivalent currents us- 
ing a calibration of shunt voltage versus current per- 
formed prior to these tests. The computer software 
continuously updated and displayed the last five val- 
ues and the average values of both shunt voltages and 
magnet currents on the screen. These values were 
also recorded as the model was unloaded to check for 
hysteresis. The drag electromagnet current was plot- 
ted against the load to determine a calibration curve 
for the model. 

To record wind-on drag measurements, the model 
was suspended at the same location where the cali- 
bration was performed. For a given Mach number, 
shunt voltages and magnet currents as well as tun- 
nel conditions ( Ps,pt , and T$) were recorded with the 
DAU and the microcomputer. A typical data point 
contains the last five readings of the voltages, cur- 
rents, and tunnel parameters. 

The drag force was obtained directly from the 
drag electromagnet current through the use of the 
calibration curve fit shown in figure 14. Buoy- 
ancy corrections were applied to this drag-force data. 
These corrections were rather large, primarily be- 
cause of the low drag forces present on the laminar- 
flow body. The drag coefficients corrected for buoy- 
ancy were reduced by 7 to 14 percent, depending on 
the tunnel speed. 

Both solid and wake blockage corrections were 
also applied to the results of this study. The drag 
coefficients corrected for blockage were reduced by 
5 percent. Further details on the data reduction and 
corrections are in the appendix. 
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Discussion of Experimental Results 

Drag-Force Calibration 

As mentioned previously, the calibration was per- 
formed with the setup shown in figure 12. Electro- 
magnet current measurements were recorded as the 
model was loaded from 5 to 242.5 gf. Current mea- 
surements were also recorded as the model was un- 
loaded to check for hysteresis. The results from this 
calibration, shown in figure 14, indicate that there 
is no significant hysteresis in these data. The data 
were fitted with a second-order polynomial. Fig- 
ure 15 shows an error plot for this drag-force cali- 
bration. These errors are the differences between the 
applied loads and the calculated loads from the poly- 
nomial equation. The root- mean-square (rms) value 
for these errors is 0.27 gf, or 0.11 percent of the full- 
scale calibration load of 242.5 gf. 

The small uncertainty of ±0.11 percent in this 
calibration may be attributable to a small amount 
of magnetic hysteresis present in the model core. 
Another possible cause is intermittent vibrations of 
the model due to slight vibrations of the optical 
position sensing system. Prior to these tests, a study 
was performed in the 13-in. MSBS comparing the 
effects of core material on calibration results. Two 
different core materials were used, mild steel and 
the same ultralow-carbon iron used with the present 
study. The ultralow-carbon iron exhibited much less 
magnetic hysteresis than the mild steel. 

Following the wind-on data runs for the laminar- 
flow body, an additional drag-force calibration was 
performed with the same loading range repeated. 
Figure 16 shows the results of this calibration com- 
pared with the results of the first calibration from fig- 
ure 14. This second calibration was also fitted with 
a second-order polynomial. The two curve fits are 
slightly displaced, with a maximum shift of 0.92 per- 
cent. These differences can be attributed to the fol- 
lowing two factors: 

1. Inaccuracies were present in positioning the 
model at the same location for both calibra- 
tions; the inaccuracies existed primarily in the 
yaw attitude and lateral position. Recall that 
the yaw attitude and lateral position were ac- 
curate to ±0.6° and ±0.1 in., respectively. 

2. Inaccuracies were present in the pitch angle of 

the calibration line used to pull on the model. 

This angle, which was also measured by the 

cathetometer, was approximately ±0.75° from 

0° in pitch. 


Drag-Force Characteristics and Repeatability 

of Data 

Wind-on drag-force measurements were taken for 
this body at Mach numbers from approximately 0.05 
to 0.20. The results, presented in figure 17, indicate 
that C£) wet decreases as Ri increases. Also shown 
in figure 17 is a repeat run for comparison; there 
is excellent repeatability of the data. The rms of 
the differences in drag force between the first and 
second runs is approximately 1.4 gf, or 0.58 percent 
of the maximum calibration load of 242.5 gf. The 
differences are larger at the higher Mach numbers 
because of vibrations exhibited by the model at these 
speeds as a result of static aerodynamic stability 
problems. Table II lists both the corrected and 
uncorrected drag coefficients for run 1. 

In reference 28, Hoerner shows the general trends 
of drag coefficient versus Reynolds number for 
streamlined bodies of revolution. (See fig. 22 in chap- 
ter VI of ref. 28.) In the Reynolds number region 
from 10 1 2 * 4 to 10 5 , Hoerner has estimated the drag coef- 
ficient assuming there is laminar boundary-layer sep- 
aration without reattachment, as in the case of flow 
over a spheroid. From R = 10 5 to 10 6 , C D>wet de- 
creases drastically because of the onset of transition 
and turbulent flow. The turbulent flow energizes the 
boundary layer and carries the flow farther around 
the body, the result being a smaller wake and thus 
a lower drag coefficient than is present with laminar 
separation. From R = 10 5 to around 10 7 , transi- 
tion is moving forward. At around R = 10 7 , there is 
fully turbulent flow over the bodies. The data from 
the present study follow some of the same trends as 
those shown in the figure in reference 28. In particu- 
lar, the data from this study follow a nearly constant 
slope for Reynolds numbers from 4 x 10 5 to 8 x 10 5 . 
This trend is very similar to the data trend shown 
by Hoerner for Reynolds numbers from 10 4 to 10 5 . 
Also, the data from this study drop off drastically 
for Reynolds numbers from 8.0 x 10 5 to 1.3 x 10 6 . 
This trend is very similar to the data trend shown 
by Hoerner for Reynolds numbers from 10 5 to 10 6 . 

Figure 18 is an oil flow photograph from the 
present study showing a laminar separation bubble 
at R L = 1.2 x 10 6 . Figure 19 is a liquid crystal flow 
visualization photograph also showing this laminar 
separation region at the same R £. A laminar separa- 
tion line and turbulent reattachment line are shown 
in both figures. At lower values of R ^ it was observed 
that the separation line was at a more forward sta- 
tion on the body. It is estimated that this separation 
line originated between x/L = 0.50 and x/L = 0.55, 
and it was observed that this separation line moved 
rearward as R L increased. It was also observed that 
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the reattachment line originated near the tail of the 
body at a lower R L . This reattachment line moved 
forward as R L increased. Thus, this separation bub- 
ble decreased in size as Rl increased. 

Based on these flow visualization observations 
and the drag coefficient results, the following deduc- 
tions are made: 

1. As Ri increased up to approximately 7 x 10 5 , 
the laminar separation line moved very slowly 
rearward, with either no reattachment or reat- 
tachment close to the tail. The increasing 
length of the attached laminar boundary layer 
would decrease the skin-friction drag. This, 
combined with the reduction in wake size due 
to the rearward moving separation line, may 
have accounted for the gradual decrease in 
Cc.wet i n this region. 

2. For Ri > 7x 10 5 , the reattachment line moved 
rapidly forward from the tail as the separation 
line moved slowly rearward. The values of 
Cd , wet subsequently decreased significantly as 
R l increased through this region. This drastic 
decrease may have been attributable to the 
decreasing size of the separation bubble or the 
decreasing size of the wake or both. 

The development and progression of the lami- 
nar separation bubble in this study was consistent 
with separation bubbles on airfoils as noted in ref- 
erence 29. In reference 29, Roberts lists four possi- 
ble flow regimes for airfoils with separation bubbles. 
These are 

1. complete laminar separation at low values of 

Rl 

2. long bubble region as Ri increases 

3. short bubble region as Rl increases further 

4. transition prior to separation as Rl increases 
further 

As noted, the first three regimes were present during 
this study. 

References 29 and 30 point out that the velocity 
and pressure distributions on a body are increasingly 
affected as the bubble size increases. In addition, 
reference 31 points out that the turbulent boundary 
layer behind a separation bubble is thicker than the 
turbulent boundary layer for the case where there is 
smooth, attached flow without a bubble. This was 
also pointed out by Holmes (private communication 
from B. J. Holmes, NASA Langley Research Center, 
Hampton, Virginia). These factors help to substan- 
tiate the deduction above which relates the drastic 
decrease in C£> )Wet to decreasing bubble size. 


Comparison With Other Data 

As mentioned previously, this shape was also 
tested at the United States Naval Academy (refs. 17 
and 18). For the present study, the desire was to 
test the model in the 13-in. MSBS up to M = 0.5 
to overlap the lower Reynolds number data from the 
Naval Academy. However, static aerodynamic stabil- 
ity problems exhibited by the model would not allow 
testing above M = 0.2. As a result, only one data 
point from this study overlaps the data from refer- 
ence 18. However, the data may be compared by 
looking at the relative drag coefficient levels shown 
in figure 20. The tabular data from reference 18 are 
used for this comparison. Note the large discrep- 
ancy in C D ^ e t at Rl = 1 x 10 6 . This discrepancy 
may be attributed to the fact that the data in ref- 
erence 18 were not corrected for cavity pressure. In 
reference 18, Hoyt states that the corrections were 
not performed; however, he gives the necessary in- 
formation and pressure measurements to correct the 
data for cavity pressure. Hoyt suggests using the fol- 
lowing equation: 

D c = D m — Pcv-Atp 

where D m and D c are the measured and corrected 
drag forces and At p is the area inside the tail piece, 
which was 45 in 2 . The data from reference 18 were 
corrected with the above equation and the results are 
shown in figure 21. The corrected data are in good 
agreement with the results from the present study. 
The difference that still exists at Rl = 1 x 10 6 may 
be due to one or more of the following: 

1. The sting at the rear of the model used at the 
Naval Academy may have acted as an exten- 
sion of the body. This would have resulted in 
a higher fineness ratio which would have re- 
duced the overall drag coefficient. 

2. The sting divided and altered the wake region 
at the tail of the model, possibly reducing the 
overall drag. 

3. The rear portion of the model used at the 
Naval Academy was slightly modified to ac- 
commodate the sting. This model did not 
have the slight flair at the tail, as was the case 
for the 13- in. MSBS model (see fig. 9). This 
slight modification may have altered the wake 
pattern. 

Drag coefficient data on this laminar-flow body 
are also presented in reference 17. However, these 
drag coefficients do not agree with those in refer- 
ence 18. The differences are not presently under- 
stood. 
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Comparison of Free and Fixed Transition 

Several runs with fixed transition were performed 
in an effort to assess the effect of fixing transition on 
the laminar separation bubble (ref. 32). Transition 
was fixed at two different locations with different grit 
sizes. Three-dimensional abrasive grain transition 
grit was used, with transition fixed at x/L = 0.17 
and 0.50. The sizes of grit used were no. 100, no. 80, 
no. 60, and no. 40, which correspond to nominal 
grit diameters of 0.006, 0.008, 0.012, and 0.020 in., 
respectively. Reference 33 was used to properly size 
the grit. Based on this technique, no. 60 and no. 80 
grits were chosen for the x/L = 0.17 station and 
no. 40 and no. 80 grits were chosen for the x/L = 0.50 
station. In addition, no. 40 and no. 100 grits were 
used at x/L = 0.17 for comparison. 

Figure 22 shows drag characteristics for the body 
with free and fixed transition. The fixed-transition 
data are for the no. 60 grit at x/L = 0.17. At the low- 
est Reynolds number, the boundary layer was thick, 
and the grit had no effect. However, as Rp increases 
to 5 x 10 5 ,C/j iWet decreases drastically. This was a 
transitional region where the boundary layer was not 
yet fully turbulent beyond the grit. However, there 
was enough energy imparted from the grit to keep the 
boundary layer attached beyond the laminar separa- 
tion line. Rearward movement of the separation line 
resulted in a lower Cp>^ et . For R L > 1 x 10 5 , the 
boundary layer was fully turbulent beyond the grit 
and the turbulent separation point moved much more 
slowly rearward with increasing R L . This resulted 
in a much slower decrease in Cp <viet with increasing 
Rl- Figure 23 shows liquid crystal flow visualization 
results for this fixed transition case. Note the tur- 
bulent separation stations. The location was farther 
back on the model for the higher Reynolds number. 
As a comparison, the laminar separation station from 
figure 19 was at x/L = 0.63. 

Figure 24 shows drag characteristics for fixed 
transition at x/L = 0.17 for all four grit sizes. As 
the grit size decreases (grit number increases), the 
significant reduction in Cp wet occurs at successively 
higher values of Rl because the smaller grit was less 
effective in fully tripping the thicker boundary layer 
at the lower Reynolds numbers. As R L increased, 
the boundary layer thinned and the smaller grit 
effectively tripped the boundary layer to turbulent 
flow. 

Figure 25 shows drag characteristics for free tran- 
sition and for transition fixed at x/L = 0.50 with 
two grit sizes. Transition was fixed at this location 
in order to retain laminar flow over a large portion 
of the body while laminar separation was avoided by 
making the boundary layer turbulent just ahead of 


the separation bubble. It was also envisioned that 
tripping the boundary layer in this region would al- 
low the turbulent boundary layer to remain attached 
farther on the body, thus resulting in a lower CJ p wet . 
The data in figure 25 show the same trends as those 
shown in figure 24. The smaller grit (no. 80) did 
not trip the boundary layer to fully turbulent flow 
until Rp = lx 10 6 . The larger grit (no. 40) ap- 
pears to have effectively tripped the boundary layer 
through the entire Reynolds number range. In fact, 
at Rp > 8 x 10 5 the boundary layer had thinned to 
the point where the no. 40 grit was too large. This 
resulted in grit drag for Reynolds numbers greater 
than 8 x 10 5 , with Cp) y, et being larger for no. 40 grit 
than for no. 80 grit. 

Comparison With and Without Fins 

A set of fins were installed in the tail of the model 
(see fig. 26) to assess their effectiveness in overcoming 
the static aerodynamic instabilities exhibited by the 
model at the higher Mach numbers. The fins did al- 
low testing at higher speeds (up to M = 0.3); thus, it 
was verified that the instabilities were aerodynamic. 
The increase in drag coefficient for the model with 
fins compared with that for the model without fins 
is shown in figure 27. The change in C D wet is rela- 
tively constant except at the higher Reynolds num- 
bers, where the increase is greater. 

Figure 28 shows drag characteristics for free and 
fixed transition for the finned model. Transition was 
fixed at x/L = 0.50 with no. 80 grit. A comparison of 
these data with fixed-transition data (same grit size) 
from figure 25 shows large differences at the higher 
Reynolds numbers. These differences may have been 
due to an altered flow pattern for the finned model 
as a result of the interaction of the fins with the thick 
turbulent boundary layer at the model tail. 

Comparison of Experimental and 
Computational Results 

SANDRAG Code 

A comparison of the drag coefficient results from 
SANDRAG with the free-transition data generally 
indicates a large discrepancy between theory and ex- 
periment. (See fig. 29.) Further insight into the com- 
putational results is given in figure 30, which shows 
a breakdown of the various drag coefficient compo- 
nents and their contributions to the overall drag co- 
efficient. For a given shape, SANDRAG breaks down 
the overall drag into four components: skin-friction 
drag, base pressure drag, pressure drag, and separa- 
tion drag. The base pressure drag is simply the com- 
puted base pressure coefficient multiplied by the ratio 


8 


of the base area to the reference area. The pressure 
drag is the drag based on the pressure distribution 
around the portion of the body that has unseparated 
flow. Note that there is not a separation drag com- 
ponent for the laminar-flow body in figure 30. At 
the lower Reynolds numbers (Rl < 7 x 10 5 ), part 
of the discrepancy in figure 29 may be attributable 
to the lack of separation drag in the computation. 
This is the drag based on the pressure in the sepa- 
rated flow region. Recall from previous discussions 
that the boundary layer may have separated com- 
pletely in this Reynolds number range. Also note 
that the experimental and computational results ap- 
proach each other as Rl increases, as shown in fig- 
ure 29. This may be attributable to the decreasing 
bubble size with increasing Rl ■ From previous dis- 
cussions, it appeared that Cx> iWe t decreased as the 
bubble size decreased because of a decreasing sepa- 
ration drag component. 

The code did predict a separation bubble on the 
body in the turbulent-boundary-layer region follow- 
ing transition. However, it did not calculate a pres- 
sure drag contribution from this bubble. (Separa- 
tion drag is only calculated when the flow does not 
reattach.) It appeared from the flow visualizations 
that the separation point was in the laminar region, 
in which case transition actually occurs on the bub- 
ble (ref. 30). Boundary-layer transition on the bub- 
ble would result in a thicker boundary layer in the 
turbulent region when compared with transition oc- 
curring prior to the bubble. This thicker bound- 
ary layer would result in a larger momentum deficit 
at the model tail and, consequently, a higher drag 
coefficient. 

The code also predicted that the bubble would 
decrease in size as Rl increased. This trend was 
also observed with the oil flow visualization. Fig- 
ure 31 shows the predictions of turbulent separation 
and reattachment points for this bubble. Note the 
differences between the locations of the experimental 
and computational bubbles by comparing the oil flow 
results from figure 18 with the computational results 
in figure 31. At R L = 1.2 x 10 6 , the oil flow sepa- 
ration line was at x/L = 0.67 and the reattachment 
line was at x/L = 0.73. The computational sepa- 
ration point was x/L — 0.835 and the reattachment 
point was x/L = 0.878. 

Nakayama and Patel Code 

As noted previously, this code calculates the skin- 
friction portion of the drag. It does not account 
for pressure drag, such as separation and base drag. 
Because of this, the predicted C/y wet values are low, 
as shown in figure 32. These predicted C jj.v/ei values 


are very close to the skin-friction drag coefficients 
from SANDRAG shown in figure 30. 

For this study, two different empirical correlations 
incorporated in the code were used to predict tran- 
sition locations. One correlation is from Granville 
(ref. 34) and the other is from van Driest and Blumer 
(ref. 35). With Granville’s criterion, this code pre- 
dicted laminar separation prior to transition, with 
a laminar separation location of x/L = 0.67 pre- 
dicted for the entire Reynolds number range. The 
code assumes immediate reattachment of the turbu- 
lent boundary layer. With van Driest and Blumer’s 
criterion, the code predicted transition at x/L = 0.66 
for the entire Reynolds number range. The code did 
not predict separation with this criterion. As shown 
in figure 32, both criteria give essentially the same 
drag coefficient results. 

A comparison is also made with some of the fixed- 
transition data by incorporating the experimental 
transition location into the code. This comparison is 
shown in figure 33, where the data are from figure 25 
(no. 40 grit at x/L = 0.50) and the computational re- 
sults are based on a transition location of x/L = 0.50. 
The predicted Cx> )We t values are close to the experi- 
mental data. The differences are greater at the higher 
Reynolds numbers because of increasing grit drag for 
the experimental data. The computational results 
are fairly accurate in this case even though the code 
calculates skin- friction drag only. If the flow was fully 
attached and if the base drag contribution was small, 
then skin-friction drag would dominate and a good 
correlation between theory and experiment would be 
expected. However, the flow is not fully attached in 
this case, as is shown in figure 34, which is a liquid 
crystal flow visualization photograph showing turbu- 
lent separation at x/L = 0.80. One possible scenario 
is that the drag for this fixed-transition case is so 
dominated by skin-friction drag that the base and 
separation drag components are very small in com- 
parison. 

Conclusions 

The primary objective of this research was met by 
substantiating the drag-force measuring capabilities 
of the 13-in. Magnetic Suspension and Balance Sys- 
tem. The drag-force calibrations were very repeat- 
able, with no significant hysteresis. The root mean 
square (rms) of the errors for the initial calibration 
was small at 0.11 percent of the full-scale calibration 
load. Wind-on drag-force measurements showed ex- 
cellent repeatability, with the rms of the differences 
between two sets of data at only 0.58 percent of the 
full-scale calibration load. 

The free-transition values of drag coefficient 
Cd , wet f° r the laminar-flow body generally decreased 
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with increasing Reynolds number R L . The drastic 
drop in Cjy ^ e t at the higher Ri values appears to 
have been directly related to the decreasing size of 
the laminar separation bubble through this Reynolds 
number region. Comparisons of these data with data 
on general streamlined bodies showed similar trends. 
Comparisons with other data (corrected for cavity 
pressure) for the same shape showed a consistent level 
and trend with Reynolds number. 

Tests with transition fixed generally showed a 
reduction in drag coefficient because of elimination of 
the laminar separation bubble. Fixing transition at 
the 50-percent model length station with no. 80 grit 
gave the lowest C value for these investigations. 
A further reduction in Cp )Wet may be obtained by 
fixing transition just ahead of the laminar separation 
line. This lower drag coefficient would result from 
a longer run of laminar flow on the body prior to 
transition. 

The model exhibited aerodynamic instabilities 
that did not allow testing above a Mach number of 
0.2. A set of fins installed in the model tail allowed 


testing up to a Mach number of 0.3; thus, it was 
verified that the instabilities were aerodynamic. As 
expected, the drag coefficients with fins were higher. 

Surface flow visualization through the use of both 
liquid crystals and pigmented oil flow verified the ex- 
istence and location of the laminar separation bubble. 
The oil flow was especially useful in highlighting the 
separation and reattachment lines. 

Comparisons with two simple computational drag 
prediction schemes were generally poor. These codes 
were developed to predict drag for simple axisymmet- 
ric shapes; they were not designed to accurately pre- 
dict drag for complicated axisymmetric shapes with 
separated flow, such as the one used in this study. 
Proper modeling of the laminar separation bubble 
and the transition location on this bubble are essen- 
tial to achieve accurate drag prediction in this case. 

NASA Langley Research Center 
Hampton, VA 23665-5225 
February 23, 1989 
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Appendix 

Data Reduction and Corrections 

Drag-Force Equation 

When a soft iron core is used, a drag-force cal- 
ibration curve fit from the 13-in. MSBS takes the 
following form: 

D = AI 2 + BI + C (Al) 

The coefficients A,B , and C change for different 
models or cores and/or different locations of the 
model or core in the electromagnetic fields. 

During a typical drag-force calibration for this 
study, the first data point was taken with a 5-gf 
weight attached to the calibration line. To get the 
drag electromagnet current value for a zero axial 
load, equation (Al) was solved with the quadratic 
formula. This formula is defined as 

-b±Vb^J ac 

2A 

The positive root J 0 is the drag electromagnet current 
value for a zero axial load. 

When the wind-on drag currents are converted to 
equivalent forces, a delta current is initially obtained 
for each data point from 

A/ = / W o — Avoff,0 (A3) 

The corrected current for each data point is then 
defined as 

/c = / 0 + A/ (A4) 

The drag force was computed with I c used in equa- 
tion (Al). This corrected current was used to com- 
pensate for possible wind-off current zero shifts. 

Buoyancy Corrections 

Three-dimensional buoyancy corrections were ap- 
plied to the drag- force data with the following equa- 
tion from reference 36: 

AD B = - < A5 > 

The A term for a body of revolution is defined as 

* = 4 jptf* (AO 


This constant term was calculated from the po- 
tential solution results from the computer code 
SANDRAG. The value of A for the laminar-flow body 
was 2.6994. For the term, tunnel-empty pressure 
measurements were taken along the centerline of one 
of the vertical test section walls. This pressure gra- 
dient varied with tunnel speed; therefore, a curve fit 
of ^ as a function of M was used when these cor- 
rections were applied to the data. 

Blockage Corrections 

Both solid and wake blockage corrections (ref. 37) 
were applied to the results of this study. The block- 
age factor is defined as 

ej> = € s i + e w (A7) 

The solid blockage term ey for three-dimensional 
models in rectangular test sections is defined as 



The r term is defined as 

1 oo oo / .. \ 3/2 

£ £ (sot <A9) 

m=—oo n=— oo N 

except for (m,n) = (0,0). This equation gives a 
value of t for a rectangular test section of width b 
and height h. For this investigation, the test sec- 
tion shape was a modified octagon which can be 
considered a rectangle with fillets. In reference 38, 
Batchelor showed that the value of r for a 7- by 
9-ft rectangular test section with fillets (similar to 
the shape used with this investigation) is very close 
to the value of r for the basic rectangle. Therefore, 
for this study the value of r was assumed to be the 
value of the basic rectangle. The basic rectangle is 
12.56 by 10.69 in. and r is 0.818. 

The wake blockage term e w for three-dimensional 
models in rectangular test sections is defined as 

<Al0) 

For this study, equations (A8) and (A10) were used 
with bh replaced by the cross-sectional area of the 
test section. This modification is suggested in ref- 
erence 37 as an approximation for octagonal test 
sections. 
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Table I. Model Coordinates for Laminar-Flow Body of Revolution 


x, in. 

r, in. 

0.0000 

0.0000 

.0010 

.0163 

.0020 

.0231 

.0040 

.0327 

.0060 

.0400 

.0080 

.0462 

.0100 

.0517 

.0150 

.0634 

.0200 

.0732 

.0250 

.0819 

.0300 

.0898 

.0400 

.1037 

.0500 

.1161 

.0600 

.1273 

.0700 

.1377 

.0800 

.1474 

.0900 

.1565 

.1000 

.1651 

.2000 

.2359 

.3000 

.2917 

.4000 

.3399 

.5000 

.3832 

.6000 

.4232 

.7000 

.4605 

.8000 

.4958 

.9000 

.5293 

1.0000 

.5614 

1.2000 

.6218 

1.4000 

.6781 

1.6000 

.7308 

1.8000 

.7805 

2.0000 

.8274 

2.2000 

.8717 

2.4000 

.9136 

2.6000 

.9532 

2.8000 

.9907 

3.0000 

1.0259 

3.2000 

1.0592 

3.4000 

1.0904 

3.6000 

1.1196 

3.8000 

1.1468 

4.0000 

1.1722 

4.2000 

1.1956 

4.4000 

1.2172 

4.6000 

1.2369 

4.8000 

1.2547 

5.0000 

1.2707 

5.2000 

1.2849 

5.4000 

1.2972 

5.6000 

1.3077 


x, in. 

r, in. 

5.8000 

1.3164 

6.0000 

1.3243 

6.2000 

1.3284 

6.4000 

1.3317 

6.6000 

1.3332 

6.6660 

1.3333 

6.8000 

1.3329 

7.0000 

1.3297 

7.2000 

1.3225 

7.4000 

1.3100 

7.6000 

1.2910 

7.8000 

1.2645 

8.0000 

1.2300 

8.2000 

1.1868 

8.4000 

1.1348 

8.5000 

1.1055 

8.6000 

1.0740 

8.7000 

1.0404 

8.8000 

1.0048 

8.9000 

.9673 

9.0000 

.9280 

9.1000 

.8872 

9.2000 

.8449 

9.3000 

.8015 

9.4000 

.7571 

9.4500 

.7346 

9.5000 

.7120 

9.5500 

.6894 

9.6000 

.6667 

9.6500 

.6440 

9.7000 

.6215 

9.7500 

.5994 

9.8000 

.5778 

9.8500 

.5569 

9.9000 

.5368 

10.0000 

.4994 

10.1000 

.4663 

10.2000 

.4378 

10.3000 

.4141 

10.4000 

.3953 

10.5000 

.3811 

10.6000 

.3713 

10.7000 

.3654 

10.8000 

.3629 

11.0000 

.3659 

11.2000 

.3752 

11.4000 

.3861 

11.6000 

.3947 

11.8000 

.3992 

12.0000 

.4000 


Table II. Drag Coefficients for Laminar-Flow Body of Revolution (Run 1) 


Rl 


Uncorrected 


3.645 x 10 5 


1.333 x 10 


-2 


4.376 

5.047 

5.608 

6.894 

7.290 

7.972 

8.482 

8.964 

9.406 

10.071 

10.775 

11.159 

11.538 

12.199 

12.575 

12.979 


1.351 

1.297 

1.259 

1.181 

1.152 

1.094 

1.044 

.972 

.900 

.783 

.698 

.647 

.623 

.583 

.548 

.562 


W}, wet 


Corrected 

1.179 x 10 -2 
1.201 
1.155 
1.121 
1.050 
1.023 
.969 
.922 
.854 
.786 
.675 
.594 
.546 
.523 
.485 
.451 
.465 
.451 


• Afterbody drag 



• Trim and 
tail loads 



't 


• Yaw damping 




a, deg 


Figure 1. Examples of model support interference problems. 
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(b) Photograph. 

Figure 2. Arrangement of 13-in. MSBS. 



L-86-1 1,940 

Figure 3. View from control room of 13-in. MSBS. 



L-86-1 1,936 

Figure 4. Laminar-flow body suspended in 13-in. MSBS. 
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Figure 5. View from contraction section of wind tunnel of laminar-flow body suspended in 13-in. MSBS. 
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Figure 8. Laminar-flow model used in 13-in. MSBS. 


L-87-11,132 
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Figure 9. Tail modification of Navy model from design shape. 
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L-87-8070 


Figure 11. Optical cathetometer on tracks. 
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mum 



L-87-8941 

Figure 13. Control room used for data acquisition. 



D, gram— force 

Figure 14. Drag-force calibration for laminar-flow body. 
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Drag— force errors, gram— force 


1.0 

.8 

.6 

.4 

.2 

0 

-.2 

-.4 

-.6 

-.8 

- 1.0 


□ 


[O 


o 


O Increasing load 
□ Decreasing load 


□ 


□ 


0 o 


□ □ 


TT 


□ 


O 


tr 


□ 


o 


o 


o 


o 


-B- 


0 50 100 150 200 

Applied drag force, gram— force 


O 


O 


250 


Figure 15. Drag-force errors from calibration of laminar-flow body. 



Figure 16. Repeatability of drag-force calibration for laminar-flow body. 



.3 .4 .5 .6 .7 .8 .9 1 2 x 10 6 

R L 

Figure 17. Drag characteristics of laminar-flow body. 



Laminar separation (x/L = 0.67) 


Turbulent reattachment (x/L = 0.73) 


L-87- 10,200 

Figure 18. Oil flow visualization for free transition. Ri = 1.2 x 10 6 . 
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Turbulent reattachment (x/L = 0.74) 



Laminar separation (x/L = 0.63) 


L-89-39 

Figure 19. Liquid crystal flow visualization for free transition. Photograph artificially enhanced because of 
poor contrast on original negative. R i = 1.2 x 10 6 . 



Rl 

Figure 20. Measured drag characteristics and uncorrected data from reference 18. 
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Figure 21. Measured drag characteristics and data from reference 18 corrected for cavity pressure. 
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Figure 22. Drag characteristics of laminar-flow body with free transition and transition fixed at 
17-percent station with no. 60 grit. 



(a) R l = 7.6 x 10 5 . 



(b) R l = 1.2 x 10 6 . 


L-89-40 


Figure 23. Liquid crystal flow visualization for fixed transition at x/L = 0.17 using no. 60 grit. Photograph 
artificially enhanced because of poor contrast on original negative. 
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.02 


Grit no. 



Figure 24. Drag 
sizes. 


characteristics of laminar-flow body with fixed transition at x/L = 0.17 for four different grit 
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Figure 25. Drag characteristics of laminar-flow body with free transition and transition fixed at 
50-percent station with both no. 80 and no. 40 grits. 
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Figure 26. Laminar-flow model with fins. Dimensions are in inches unless otherwise noted. 
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Figure 27. Drag characteristics of laminar-flow body with and without fins. 
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Figure 28. Drag characteristics of laminar-flow body with fins for free transition and transition fixed at 
x/L = 0.50 with no. 80 grit. 


31 


.02 


□ SANDRAG results 
O 13— in. MSBS data 
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Figure 29. Drag characteristics of laminar-flow body computed with SANDRAG and measured during free 
transition. 



Rl 

Figure 30. Breakdown of SANDRAG drag components. 
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Figure 31. Separation and reattachment points from SANDRAG and from oil flow results. 
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Figure 32. Drag characteristics of laminar-flow body computed with Nakayama and Patel code and measured 
during free transition. 
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Figure 33. Drag characteristics of laminar-flow body computed with Nakayama and Patel code based on 
experimental transition location and measured for fixed transition at x/L = 0.50 with no. 40 grit. 



Figure 34. Liquid crystal flow visualization for transition fixed at x/L = 0.50 with no. 40 grit. R l = 1.1 x 10 6 . 
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